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In this paper, we investigate the fermions Hawking radiation and quasinormal modes in infra-red 
modified Hofava-Lifshitz gravity under tunneling and perturbation perspectives. Firstly, through 
the fermions tunneling in IR modified Hofava-Lifshitz gravity, we obtain the Hawking radiation 
emission rate, tunneling temperature and entropy for the Kehagias-Sfetsos black hole. It is found 
that the results of fermions tunneling are in consistence with the thermodynamics results obtained 
by calculating surface gravity. Secondly, we numerically calculate the lowing quasinormal modes 
frequencies of fermions perturbations by using WKB formulas including the third orders and the 
sixth orders approximations simultaneously. It turns out that the actual frequency of fermions 
perturbation is larger than in the Schwarzschild case, and the damping rate is smaller than for the 
pure Schwarzschild. The resluts of fermions perturbation suggest the quasinormal modes could be 
lived more longer in Hofava-Lifshitz gravity. 

PACS numbers: 04.70.Dy, 04.62. +v, 03.65.Sq 

Keywords: Horava-Lifshitz gravity; fermion tunneling; quasinormal modes 



I. INTRODUCTION 



Recently, Hof ava presented a power counting renormalizable gravity theory at Lifshitz point which is called Hofava- 
Lifshitz (HL) gravity It exhibits a broken Lorentz symmetry at short distances and reduces to usual general 
relativity (GR) gravity at the large distances with particular A = 1 which controls the contribution of the extrinsic 
curvature trace. With HL gravity theory putting forth, HL gravity is intensively investigated in many aspects involving 
basic formalism Q, cosmology [j], various black hole solutions and their thermodynamics 0-1? | and so on. 

In the subsequent developments of the HL gravity, people are trying to find the influence of matter fields for HL 
gravity. Which analogue of the matter energy-momentum tensor could be used to act gravitational source as GR? The 
jioneering works involve the geodesic analysis by various methods including the optical limit of a scalar field theory 



Relations 
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, super Hamiltonian formalism [9] , foliation preserving diffeomorphisms [lol | , Lorentz- violating Modified Dispersion 



an so on. The optical limit presented by Capasso and Polychronakos [8] could offer a deformed geodesic 
equation by the generalized Klein-Gordon action. It is shown that the particles maybe does not move along geodesic. 
Deviations from geodesic motion appear both in fiat and Schwarzschild-like spacetimes. Similar result deviations from 
GR also is found in [s-Hl with various above-mentioned approaches. 
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Considering vanishing cosmological constant (Aw), Kehagias and Sfetsos proposes an asymptotically flat black 
hole solution by introducing the addition term proportional to the Ricci scalar of three geometry fi^R^^\ which indi- 
cates the Minkowski vacuum and modified GR at infra-red (IR) modification. Then after that, many people have de- 
voted to its phenomenolo gy i nvolving strong field gravitational lensing[l^[3; scalar field quasi-normal modes [l2l.[l^. 
timelike geodesic motion [isl , thin accretion disk 16 [ and observations constraints as well as its thermodynamics 
analysis |18l420l |. which presents the black hole entropy S = A/4 + 7r/aln(j4/4) via the first law of thermodynamics 
where a is Hofava parameter. This entropy could be treated as Generalized uncertainty principle quantum correction 
entropy [3] or casting entropy If fermions are tunneling in this kind HL gravity, might it keep this logarithm 
entropy? Motivated by this, we investigate fermions tunneling in section III. 

On the other hand, as we know that through some additional fields (e.g. scalar or fermions fields), the black hole 
suffers a damping oscillation phase which is named as "quasi-normal models" (QNMs) or "quasi-normal ringing" . As 
a results, the normal model oscillation is replaced by a complex frequencies which encodes the black hole's important 
information such as mass, charge, momentum and the dimensions of spacetime. The real part of complex frequency 
represents the actual frequency and the imaginary part of its represents the damping of the oscillation. It is believed 
that these QN frequencies could be detected by (LIGO, VIRGO, TAMT, GEO600) in future. So according to this 
observable QNMs, some people have used massless scalar field to obtain the important lowing QN frequencies which 
could live longer and be detected easily Q] • It is interesting that the QNMs of massless scalar field are longer lived 
and have larger real oscillation frequency in Hofava-Lifshitz gravity than in GR. Else, as one kind of basic particles, 
the fermions could offer us many important information. Motivated by the situations above, we will evaluate the 
QNMs for the massless fermions perturbations in section IV. 

As we have already known that, due to the different kinetic terms, Hofava-Lifshitz gravity with A 7^ 1 is different 
significantly from General Relativity. Hence, the HL gravity with A 7^ 1 has been extensively studied in the literatures. 
The relevant works are mainly concentrated on the basic problems such as how to find various exact black hole 
solutions, the application of cosmology, the constraints of various fundamental parameters and so on. Despite more 
attention has been paid to the properties of HL black holes, there are a few works referring to the fermions analysis, 
especially to the A 7^ 1 case. So far as we know, the works relevant the fermion analysis focus mainly on the IR 
modified HL gravity includin g D irac perturbations 



on. In Wang and Gui's work 
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and fermion tunnelling for z = 4 black holes 



231 and so 



21j . the quasinormal frequencies of massless Dirac field perturbation are evaluated by 



22| , the evolution of Dirac perturbations is also 



23|, the 



third-order WKB approximation. In Varghese and Kuriakose's work 

investigated by using time domain integration and third-order WKB methods. In Chen, Yang and Zu's work 
fermion tunnelling is investigated in the background of (3 + 1) dimensions and (4-1-1) dimensions black holes in the 
z = 4 HL gravity. 

This paper is organized as follows. In section II, we present the Kehagias and Sfetsos black hole solutions. In section 
HI, we calculate the fermions tunneling emission rate of Hawking radiation, tunneling temperature and entropy. In 
section IV, we use the third and the sixth orders WKB formulas to numerically calculate frequencies simultaneously. 
Section V is the conclusions. We adopt the signature (— , +, +, +) and put h, c, and G equal to unity. 
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II. AN ASYMPTOTICALLY FLAT INFRA-RED MODIFIED BLACK HOLE SOLUTION IN 

DEFORMED HOrAVA-LIFSHITZ GRAVITY 

In this section, we review briefly the KS black hole solutions under the limit of K\y — > with running constant 
A = 1 in the IR critical point z = 1. The space geometric is parameterized with Arnowitt-Deser-Misner (ADM) 
formalism, 

ds^ = -N^dt^ + g,j (dx' + N'dt) [dx' + W dt) . (1) 
The action for the fields of HL theory is 



where the second fundamental form, extrinsic curvature Kij, and the Cotton tensor C^^ are given as follows, 

2N \dt" ' 



= ^ (^9^J - V.iV, - V,7V,; ) , (3) 



(4) 



Here, n, A, a, /i and Aw are the constant parameters. The last term of metric Eq.Q represents a soft violation of 
the detailed balance condition. Comparing the HL gravity action with that of GR gravity, we can obtain the speed 
of light c, the Newton's constant G and the cosmological constant A 

In the limit of Aw — > 0, we can obtain a deformed action as follows, 

S ^ [ dtLd^xiCo + Ci) , (6) 

(7) 



Co - ^nS^^{K,jK'^ -XK^)^ 



A = V^Nr^Jl^n^ - ^ ( a, - l^R.,) ( - f^ir^ ]^,^n}. (8) 



32 (1 - 3A) 2a4 y 'J 2 V V 2 

For the particular case of A = 1 with a = a spherically symmetric black hole solution is presented by 

Kehagias and Sfetsos which also is corresponding to an asymptotically flat space, 

ds^ = ~f(r)dt^ + — - + (dO^ + sin^ 0dip^) . (9) 
/('') 

The lapse function is 



/(r) = l + ar^ - r (aV3 +4aM), (10) 

where the parameter M is an integration constant related with the mass of black hole. Using the null hypersurface 
condition, one can find there are two horizons, inner r_ and outer event horizon r+ in this space. 



^± = M|1±,/1-^ . (11) 
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Thermodynamic quantities including mass Mks, temperature TkSi a-nd heat capacity Cks and entropy Sks presented 
in Refs.|18| are hsted as, 

1 + 2ar±2 



Mks = 
Tks = 
Cks — 



Aar± 
2ar\ - 1 
SiTr+{ar\ + 1)' 



27r 

a 



{ar\ + lf{2ar\ - 1) 



A 



Sks ^ + - In — 



2a'^r'\ — 5ar'i_ — 1 



A 



(12) 
(13) 
(14) 
(15) 



with horizon area A = Airr'^. Under the Hmit of a 
Sbh — A/ A for Schwarzschild black hole. 



-co, the entropy reduces to Bekenstein-Hawking entropy 



III. FERMIONS TUNNELING OF IR MODIFIED HOfAVA-LIFSHITZ GRAVITY 



In this section, we investigate the Hawking radiation of Kehagias and Sfetsos black hole in IR modified Horava- 
Lifshitz gravity with fermion tunneling. The tunneling probability, temperature and entropy are expected to be 
obtained. The Dirac equation in the KS black hole spacetime can be written as 



* = 0, 



where m is the mass of fermions. is the inverse of the tetrad e° defined by black hole metric g^jy = ?7ahej^e^ 
Minkowski metric rjat — diag(— 1, 1, 1, 1). 7" is the Dirac matrix and is the spin connection given by 

where the covariant derivative of Cbi^ is given by Christoffel symbols F^^ as, 



(16) 
^ with 



We choose following 7 matrix. 




-ia' 
2 







where cr' is Pauli sigma matrix. 




In the presentation of tx*, the spin up wave function is written as, 

ip^{t,r,e,(j)) = exp 

fA{t,r,0,<P)\ 






exp 



,0 -1 



-/t(i,r,0,, 



(17) 



(18) 



(19) 



(20) 



(21) 
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where denotes the eigenvector of spin up state with eigenvalue +1, and ^j, denotes the eigenvector of spin down 
state with eigenvalue —1. is the action of radiation particles with spin up. According to e^e^^ — (5^, we have 



e^ = diag(^, ^/f, 

\ y/j r r sm 

Submitting Eg. ip^ into Dirac Ea. ([TB|) . the frame should satisfy following relation 
Simplifying above Eq. ([23|) . we can get 



^// r rsmb' n ' 



If we neglect the small quantity F^, Ea. ([M)) could be simplified as, 



23 

r rsmo' n ' 



On the benefit of 7 metrics Eq. ((20|) . we have 

/Aiexp^a,Jt\ 


' n 1 «-ft o r 

--^nexp^^t/^ 

v y 





7^a0*t 



Submitting Eqs. 



\-iAj. exp ^9e/t/ 
and (P7)) into Dirac Eg . ([25)) . we can get 
/ Adtl-t \ ( Bdrif \ / \ 



/ Siexp^a,/t \ 


-H exp '-^drl^ 

V / 

/ \ 





1 

V7 





V y 



V7 





-Adrif 

V y 







1 



rsin^ 



/ \ 

Bd^I^ 




fA\ 


voy 



This equation could be reduced to four components of {t,r,0, 4>) as, 

A 



-B 



dtit + B/Jp)drlf + mA 



dtit - A^/f{rjdrIt + mB 



— ( idel^ 



1 



sin 9 



0, 
0, 

0, 
0. 



Consider the symmetry of the spacetime, we adopt the action below as, 

= --wi + WW + 6(61,0). 



(22) 



(23) 



(24) 



(25) 



(26) 



(27) 



(28) 



(29) 

(30) 

(31) 
(32) 

(33) 
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Submitting Eq.§^ into Eqs.(l29l), dSO]), ISI]), dSS]), we can get 

A 



B 



^JJ + B^/f{r)^rW + mA = 0, 



:uj - A^/J{r)^rW + mB = 0, 



B [^dgQ 
A ( idgQ 



1 



sin 9 
1 

sin 9 



0. 



(34) 
(35) 
(36) 
(37) 



Because the contribution of 9 on outgoing particle is equal with that of 8 on incoming particles, Eas. ([55)) and ([57)1 
do absolutely nothing that are useful to the calculation of the tunneling probability, We only need consider the action 
of the radial direction, i.e. Eas. (|M|) and (|55|). whose solvability condition is the determinant of the coefficients of A 



and B is zero. Namely, 



By direct integration of determinant Eq. 



, W could be obtained as. 



= 0. 



w±(r) = ±y" 



y^uj^ - ■m?f{r) 
fir) 



dr. 



(38) 



(39) 



Using the condition /(r) — > near horizon r+, the numerator of integrated fraction in Eq. p9|) is reduced to 
^ up' — m'^f(r) — > io. Hence, the terms contained mass (^ fn"^ f) has nothing to do with the tunneling probability. 
So, yV±[r) is applicable to the whole fermions, no matter massive or massless particles. Adopting the contour 
integration, we can get 

u 



W±{r) = ±in—r 



(40) 



where " + " denotes outgoing fermions, 
respect to r, 



denotes incoming ones, ' means the first-order derivation of /(r) with 



/'(.+ ) 2a.- 



2a^r^ + 2aM 



(41) 



dr y/ah^^+IoMr 
It is well known that the tunneling probability could be related to the imaginary part of the action. Thus, the 
tunneling probability of the emission fermion is written as followings, 

P{emission) exp(— 2/to/-|-_|_) exp(— 2/TOyV4.) 



Piahsorption) exp(— 2/m/-|-_) exp(— 2/r7iyV_) 
Submitting >V±(r) into above Ea. p2)) . we can obtain 



exp 



47rcj 

TW) 



ar^ + AM + (ar^ + Ma) ^1 + AM j Oir\ 



(42) 



(43) 



AaMr\ - 2a^r\M - a^M^ 

According to the usual relation between inverse temperature /3 and tunneling probability, F = exp(— /3a;), we can get 
the fermions tunneling temperature, 



fermion 



fir+) 
An 



1 

2^ 



a'^r^ + aM 



a^r\ + AaMr+ 



(44) 



If we choose the mass function definded by Eq. (IT2]) in Ref.jlSj, we can get 

_ 2arl - 1 

Tfermion — 7Z ; : (45) 

87rr+(l + ar^) 

which is just the Hawking temperature of the black hole in IR deformed Horava-Lifshitz gravity As a thermo- 
dynamical system, the first law dM = TdS gives black hole entropy, 

S = I T-^dM + So = [ dr+ (L^)+So. (46) 



Tdr+ 



n[rl + -\nrl] +S„, (47) 



where we adopt Mks Eq. (jl2p . If we adopt Sq = vrlnTr/a, the final logarithmic entropy obtained through fermions 
tunneling is 



S^^ + lln^. (48) 
4 a 4 ^ ' 



Based on the surface gravity defined by 



1 df 



2ar\ - 1 

^-a-v (49) 



4r+(l + ar\) 



the thermodynamic temperature Eq. (jl3p and the thermodynamic entropy Eq. ljlSp are obtained in previous researches 
It is interesting that the resluts Eqs. ([l3|) and (|15|) based on surface gravity are agreement with fermions tunneling 
results Eqs.dini) and pS)) . 



IV. FERMIONS PERTURBATIONS OF IR MODIFIED HOfAVA-LIFSHITZ GRAVITY 

In this section, we evaluate the quasinormal modes of fermions perturbation by using the third-order and sixth- 
order WKB formulas, simultaneity. In order to get the quasinormal frequencies, we should proceed from the Dirac 
Ea. (IT6)) . According to the relation of = Vabe"^, we can have e^^ = diag(— V7, l/'\/7j'';''sin^). Then, based on 
Bbi^-fj. = d^dbv — r^i/CfiQ,, the nonzero covariant derivative could be listed as, 

eoi;0 = , eio;o = ^2i;2 = -1, ei2;2 = y jr, 

e23;3 = r sin 6* cos 6*, ci^.^ = y/Jr svc? 0, 631^3 = - sin 6*, e32;3 = cos 6*. (50) 
Submitting above Eqs. ([5D|) into Eq. pT]) . we can get the spin connections as, 

rt = -J/V7'; r, = 0; Te = -^i^i^:. = -i (sin 07/7173 + cos 07 V) ■ (51) 

Considering the symmetry of frame Ea.(P^. the Dirac Eg. dTSl) could be rewritten as, 

7°e* {dt + Tt) + i^el (9, + T,,) + 7^6^ {dg + Te) + i^e^ {d^ + T^)] $ = 0, (52) 

where we adopt the massless Dirac field to simplify the perturbation problem. 

Based on the spin connections Eg. ([51]) and the anticommutation relation of 7 metrics, Eq. ((52|) could be reduced to 
a simple form as following, 

r sm a 0(p 



7° 5$ n if d i\^-f^fd 1 \^ 
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FIG. 1: the potentials Vi{r,k) (solid line) versus radial coordinate r with k = 1,2,3. Meanwhile, we also draw Schwarzschild 
case with dotted line for comparison. 



where r, 9,4') = f ^^^{r)'^{t, r, 9, (j)). Then, we could adopt an ansatz as foUowings, 

/ F{r) \ 



r\J sin 9 



Fir) 
iG{r) 
\iG{r)) 



(54) 



Submitting the ansatz Ea. ([M)) into Eg. ([5^ . we can get three equations: one equation of VL refers to variables (6*, ( 
and two equations of G{r) and F{r) refer to variable r, which are listed as, 

„3 



-.F(.) + ^^-^G(.) =0, 
or* r 



.G(0 + ^ + ^F(.) =0, 
ctr* r 



(55) 
(56) 
(57) 



where fc = — Z or Z + 1 and the coordinate transformation dr — f{r)dr^ is adopted. Eliminating F{r) (or G{r)) in 
Eas. ((56t and ((57|) . we can obtain two 2th order differential equations of G(r) (or F(r)), 

fG 
drl 



+ {uj^~Vi)F = 0, 



iu^-V2)G = 0, 



(58) 
(59) 
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where Vi and V2 are supersymmetric partners with same spectra, 

V. = ^^(nv/7M + ^^-/m) k = i + (60) 
V. = ^^(nVTM^^ + Zw) k^-i. (61) 

In the following words, we use the Eq. (j58p contained potential Vi to evaluate the quasinormal mode frequencies of 
the massless Dirac field by the third orders and sixth orders WKB approximation. Here, Vi{r, k) is plotted in Fig[T] 
which illustrates clearly that with bigger Vi is higher than that of Schwarzschild case (dotted lines). Moreover, 
the gap between KS and Schwarzschild increases greatly as increasing k, in particular near the maximum points. So 
we can expect the QNMs of fermions perturbations could be lived more longer and the actual frequencies increase 
because there are more lower potential for IR modified Horava-Lifshitz gravity. 

According to the potential Vi (r, k) Eq. ((60)) , the massless Dirac quasinormal modes in the KS black hole spacetime 
satisfies the boundary conditions, 

^{x) ^ exp(±iw), X — > ±00. (62) 

where uj — Re(a;) +2lm(a;). The real part Re(a;) determines its actual oscillation frequency and the absolute value of 
imaginary part |Im(a;)| determines the damping rate. 



TABLE I: Various low-lying overtones QN frequencies for a fixed a = 0.5. 



\k\ 


n 






3th 








6th 




Schwarzscliild(6th) 


1 








200527 


-0 


071134i 





199480 


-0 


0655111 


0.182642 


- 0.094937* 


2 








419782 


-0 


070205i 





420578 


-0 


0694681 


0.380069 


- 0.096366i 




1 





393067 


-0 


213644i 





396805 


-0 


2085701 


0.355860 


- 0.297269i 


3 








635027 


-0 


070293i 





635293 


-0 


0701741 


0.574094 


- 0.096307j 




1 





617529 


-0 


211964i 





618963 


-0 


211069i 


0.557016 


- 0. 2927174 




2 





583072 


-0 


356887i 





586030 


-0 


3536961 


0.526534 


- 0.499713i 


4 








849148 


-0 


070367i 





849265 


-0 


070339i 


0.767354 


- 0.096270j 




1 





836227 


-0 


211639i 





836870 


-0 


211407i 


0.754300 


- 0.290969i 




2 





810528 


-0 


354547i 





811997 


-0 


3536791 


0.729754 


- 0.491909i 




3 





772412 


-0 


500217i 





774470 


-0 


498069i 


0.696728 


- 0.702337j 


5 





1 


062830 


-0 


0704101 


1 


062890 


-0 


070401i 


0.960293 


- 0.096254j 




1 


1 


052590 


-0 


211548i 


1 


052930 


-0 


211469i 


0.949759 


- 0.290149i 




2 


1 


032170 


-0 


3536631 


1 


032970 


-0 


3533521 


0.929491 


- 0.488114i 




3 


1 


001690 


-0 


4974381 


1 


002930 


~0 


496649i 


0.901073 


- 0. 692514* 




4 





961423 


-0 


643574i 





962684 


-0 


6420411 


0.866730 


- 0.905116i 



In the various methods to get the frequcies of QNMs, the WKB numerical formulas are convenient to give accurate 
frequencies values for the longer lived quasinormal models. This method is originally shown by Schutz et al 24] and is 
later developed to the third order by Iyer et al 



25 



26[. At a later time, WKB approximation of QNMs is explanded 



to the sixth order by Konoplya [27]. Then after that, this method is extensively used in various spacetimes [28[. In 
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FIG. 2: Massless Dirac quasinormal mode frequencies. For the convenience of comparison, double results including KS solution 
(solid triangle) and Schwarzschild solution (solid square) are given simultaneously. 



this paper, we numerically calculate the lowing modes frequencies through the sixth order WKB formula which has 
the form [2l\ as following, 



1^0 



- Aa - A3 - A4 - As - Ae = n + 1/2, 



(63) 



where Q is a "reverse potential" given hy Q = u —V. Qq denotes the i-th derivative of Q at its maximum point with 
respect to the "tortoise coordinate" r*. The results of the third orders could also be obtained by Ea.(|551) without A4, 
A5 and Ag. Considering WKB approximation fails to calculate the higher order modes, we only evaluate low-lying 
QNM modes (n < k) by various overtones n. The correctional terms of A2 and A3 are given in Refs. 



25 



261. The 



correctional terms of A4, A5 and Ag are given in Ref427]. It turns out that WKB series shows well convergence in 
all sixth orders for Dirac field, which is similar to the scalar field case QjQl- ^his paper, we analyse the effect of 
parameter Horava-Lifshitz gravities on QNM modes through two kinds of a: one is fixed and another is changed. 

For the fixed a, we only consider the first five low-lying modes k = 1, 2, 3, 4, 5 with < n < k. The results are 
listed in Table I. The quasinormal mode frequencies for positive k are plotted in Figl2] which illustrates the real part 
Reuj decreases with increasing mode number n for the given angular momentum number k. Else, the absolute value 
of imaginary part |/maj| increases as bigger n which indicates higher modes decay faster than the low- lying ones. 
Comparing with Schwarzschild results (solid quare points), the Reu of a is larger than Schwarzschild limit, while the 
damping rate j/mwl is smaller than pure Schwarzschild case. 

For the changed a, we treat 1/2q; changed in [0, 1] as a whole. Three important kinds low-lying modes: (k — l,n — 
0), {k = 2,n — 0) and (fc = 3, n = 0) are listed in Table II, III, and IV, respectively. According these three tables, we 
plot the real part Reuj and the imaginary part Imuj of the third and the sixth order results in Figl3] Here, it should 
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TABLE II: The QN frequencies of null overtones modes (n = 0, fc = 1). 



l/2a 






Uh 





0.176452 


- 0.100109Z 


0.182642 


- 0.094937i 


0.1 


0.179189 


- 0.098175i 


0.184899 


- 0.092152i 


0.2 


0.181871 


- 0.096140i 


0.187179 


- 0.089352i 


0.4 


0.187091 


- 0.091667i 


n 1 m /I /I o 
U. 191440 


— O.Uo39Uu 


0.5 


0.189632 


- 0.089156't 


0.193342 


- 0.081265i 


0.6 


0.192118 


- 0.086395i 


0.195055 


- 0.078615i 


0.8 


0.196809 


- 0.079810i 


0.197915 


- 0.072858i 


1 


0.200527 


- 0.071 134i 


0.199480 


- 0.06551H 



TABLE III: The QN frequencies of null overtones modes (n = 0, ft = 2). 



l/2a 








Uh 





0.378627 


- 0.0965424i 


0.380069 


- 0.096366i 


0.1 


0.381664 


- 0.0949046i 


0.383040 


- 0.094794i 


0.2 


0.384885 


- 0.0931728i 


0.386218 


- 0.093061i 


0.4 


0.391946 


- 0.0893179'/ 


0.393230 


- 0.089069i 


0.5 


0.395827 


- 0.0871212* 


0.397087 


- 0.086766i 


0.6 


0.399977 


- 0.0846759-i 


0.401200 


- 0.084209i 


0.8 


0.409203 


- 0.0786964i 


0.410286 


- 0.078038i 


1 


0.419782 


- 0.0702051i 


0.420578 


- 0.069468i 



TABLE IV: The QN frequencies of null overtones modes (n = 0, fc = 3). 



1/2q; 


3t/i 




6th 





0.573685 


- 0.096324i 


0.574094 


- 0.0963070i 


0.1 


0.578089 


- 0.094773Z 


0.578503 


- 0.0947834i 


0.2 


0.582754 


- 0.093108Z 


0.583169 


- 0.0931358i 


0.4 


0.592999 


- 0.089343i 


0.593402 


- 0.0893718i 


0.5 


0.598659 


- 0.087176i 


0.599051 


- 0.0871898i 


0.6 


0.604747 


- 0.084754i 


0.605126 


- 0.0847448i 


0.8 


0.618500 


- 0.078810i 


0.618839 


- 0.0787424i 


1 


0.635027 


- 0.070293i 


0.635293 


- 0.0701738i 



notice that the horizontal abscissa denotes the value of l/2a. We impose the interpretations on these data and draw 
conclusions from them. 

(1) The real part Reco increases as bigger l/2a and the absolute value of imaginary part |/mu;| decreases with 
increasing l/2a, which indicates these QNM could be lived longer. 

(2) The gap between the third results and the sixth order ones is visibly displayed in the imaginary part. The real 
parts of them basically have the same values, except for A; = 1 modes. In general, the average relative magnitudes of 
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FIG. 3: The real parts Reuj and the imaginary parts Imuj versus radial r under the third orders (solid triangle) and the sixth 
orders (solid pentagram) WKB approximations. 



the gaps are approximately given as, 



&th 



Im(a;) -^fi ini(^) 



3t'»Ini(cj) 
6*''Re(a;) -^fi Re(a;) 



10%, 



0. 



(64) 
(65) 



3*''Re(a;) 

Hence, the orders of WKB approxmations have the tremendous bearing on the damping rate, more than on the actual 
frequency. 

Moreover, when horizontal abscissa of Fig 13] approaches Schwarzschild case, namely 1/2q; — > (or a — > +oo), 
the real frequencies Reu decreases and the damping rate |/r7iw| increases. In another words, the Horava-Lifshitz 
gravities have the longer lived and more bigger actual frequency than that of usual Schwarzschild case. This specific 
phenomenon also is observed in the massless scalar field perturbation 

HQ 



V. CONCLUSION 

In this paper, we have investigated fermions tunneling and perturbation in the IR modified Horava-Lifshitz gravity. 
We summarize what has been achieved. 

(1) For the fermions Hawking radiation, we consider the symmetrical characteristic of spacetime and adopt the 
action with the form of Eq. (j33p . Through the decomposition of Dirac Eg. (1161) . we can get the imaginary part of 
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fermions action which could help us obtain the tunneling probability according to Eq. (|42)) . Then, on the benefit of 
r = exp(— /3a;), the tunneling temperature could be given as Eg. (j44| . So, according to the first law dAI — TdS, we 
have a tunneling entropy Ea. (l48p naturally. It is interesting that the tunneling Hawking temperature and tunneling 
entropy are agreement with that obtained by calculating surface gravity. 

(2) For the fermions perturbation, we obtain lowly damped quasinormal modes by using the sixth orders WKB 
approximations, as well as the third orders formulas. In order to get a detail analysis on these obtained quasinormal 
frequencies, we adopt two kinds of methods: one is to fix the Hofavaparameter a (fixed a) and another is to change 
a in the range [0, +oo] (varied a). 

For the fixed a case, the results turns out that the impact of the Horava-Lifshitz gravity on quasinormal frequencies 
is quite. This is profoundly manifested in the following ways: the actual frequencies becomes bigger and the damping 
rate becomes more slower which indicates these lowing modes could be lived longer than that of usual Schwarzschild. 
This fact also could be explained by the perturbation potential Vi Eq. (|60| illustrated in FiglJl which shows the 
potential contained Horava-Lifshitz gravity (solid lines) is lower than that of Schwarzschild (dotted lines). 

For the varied a case, we have calculated numerically three kinds of important lowing modes (fc — l,n — 0), 
{k — 2,n = 0) and (fc = 3,n = 1) by through the third and sixth orders WKB approximations. The result listed in 
Tables II, HI and IV show three facts as follows, (i) With bigger parameter l/2a, the real part of frequencies increases 
and the absolute value of imaginary part decreases which is illustrated by FiglSl In other words, if parameter l/2a 
becomes larger, the actual frequency of QNMs will be larger with more longer damping rate, (ii) The real part is 
not sensitive to the third or the sixth orders WKB approximates. This fact also is illustrated in the left sub-plotted 
curves in FiglSj except for some small l/2ck modes of fc = 1. (iii) Against the real part, the imaginary part is sensitive 
to our WKB approximates methods. The gap of the third and the sixth orders results is unchanged basically. This 
fact also is illustrated in the right sub-plotted curves in FiglH In a word, if these specific information could be tested 
by LIGO, VIRGO, TAMT, GEO600, it wiU support Horava-Lifshitz gravity forcefully and energetically. 
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